Flow of foam past an elliptical obstacle 
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To investigate the link between discrete, small-scale and continuous, large scale mechanical prop- 
erties of a foam, we observe its two-dimensional flow in a channel, around an elliptical obstacle. 
We measure the drag, lift and torque acting on the ellipse versus the angle between its major axis 
and the flow direction. The drag increases with the spanwise dimension, in marked contrast with 
a square obstacle. The lift passes through a smooth extremum at an angle close to, but smaller 
than 45°. The torque peaks at a significantly smaller angle, 26°. No existing model can reproduce 
the observed viscous, elastic, plastic behavior. We propose a microscopic visco-elasto-plastic model 
which agrees qualitatively with the data. 
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A foam is a model to study materials which are viscous, 
elastic, and plastic. This complex, ubiquitous behavior 
is exploited in numerous applications, such as ore sep- 
aration, drilling and extraction of oil, food or cosmetic 
industry but is not yet fully understood 0- Foam 
rheology is thus an active research area, and recent stud- 
ies provide insight to understand the interplay between 
the bubble scale and the whole foam behavior 0, 0] and 
to unify elasticity, plasticity and viscosity 0]. Here, we 
study the flow of foam around an ellipse, where the mea- 
sured lift, drag and torque show the whole complexity of 
foam rheology, which strongly constrains possible mod- 
els: simple ones do not capture the observed features. 
We propose an elastoplastic model which describes well 
the data. We discuss the generality, implications and 
limitations of this model. 

We have built a foam channel 0] to investigate a 2D 
steady flow and measure the force it exerts on an obsta- 
cle (Stokes experiment 0, 0, E3). Briefly, aim 
long, 10 cm wide tank is filled with deionized water with 
1% of commercial dish- washing liquid (Taci, Henkel). Its 
surface tension is 7 = 26.1 ± 0.2 mN m , and its kine- 
matic viscosity is v = 1.06 ± 0.04 mm s~ 2 . Several com- 
puter controlled injectors blow nitrogen in the solution 
to form a horizontal monolayer of bubbles of average 
thickness h = 3.5 mm, confined between the bulk so- 
lution and a glass top plate (quasi-2D foam) ^]J. This 
foam is monodisperse (bubble area at channel entrance: 
Aq = 16.0 ± 0.5 mm 2 ) and its fluid fraction is estimated 
to be around 7% [iH ]. 

We study here the simplest shape which symmetry 
is low enough to observe simultaneously drag, lift and 
torque: the ellipse (Fig. [JJ. The obstacle has a major 
axis 2a = 48 mm and minor axis 2b = 30 mm. It floats 
freely just below the top glass surface, without solid fric- 
tion. The upper end of an elastic fiber passes through a 
hole in the bottom of the obstacle, ensuring its free rota- 
tion, while its lower end is fixed, so that a top view of the 
obstacle displacement from its position at rest measures, 




FIG. 1: Top view of the elliptical obstacle and of the sur- 
rounding bubbles. Ellipse (a) perpendicular to the flow, (b) 
tilted in the flow, (c) parallel to the flow. The x axis is the 
direction of the flow direction and of positive drag, y is the 
direction of positive lift, and 9 (between and 90°, by symme- 
try) is the angle between x and the major axis of the ellipse. 



with a precision better than 0.1 mN, the force exerted by 
the foam on the obstacle. 

We measure the drag in the parallel orientation (9 = 
0), which is stable (see below); and in the perpendicular 
one (9 = 90°), which is unstable (but where the ellipse 
can remain for one hour, enough to perform steady flow 
measurements). The results (Fig. |2J) are very close to 
that for circles of diameters 30 and 48 mm, respectively: 
this suggests (see also Fig. 0Ji) that drag is proportional 
to the spanwise direction (along the y axis) £ of the el- 
lipse: 



= 2\J a 2 sin 2 9 + b 2 cos 2 1 



(1) 



In a steady flow (530 ml min^ 1 , i.e. a velocity of 2.5 
cm s _1 ), we start from a given initial orientation (76, 64, 
48 or 18°), let the ellipse rotate freely to its parallel stable 
orientation, and measure the angle, drag and lift (Fig. 
El- The angular velocity strongly increases in the range 
15° < 9 < 40° (with a peak at 26°) and does not depend 
on the initial orientation (insert of Fig. |3J). Moreover, 
the forces correlate to 9; we thus eliminate the time and 
plot the dependence of drag and lift with 9 (Fig. ^Ji). 
All the forces data collapse on two master curves, one for 



2 



10 



5 4 
Q 







P ^^^^ 
















200 



400 



600 



Flow rate (ml min 1 ) 



FIG. 2: Drag exerted by the flowing foam on an elliptical 
obstacle, versus the flow rate: .,0 = 0° (Fig. |TJi); ■, = 
90° (Fig. ^:). Bold lines are linear fits to the data. Data 
(open symbols) and fits (dashed lines) for circular obstacles 
(o, 30 mm; □, 48 mm diameter) from Ref. |(| are plotted for 
comparison. 



the drag and one for the lift. The drag increases roughly 
linearly with 9 except very close to 0° and 90°, where it 
is extremal by symmetry (it equals 4.5 mN for 0° and 
8.8 mN for 90°). As suggested before, the experimental 
angular dependence of the drag is close to the one of 
the spanwise dimension, despite small discrepancies for 
angles close to 0° and 90°. The lift vanishes at 0° and 
90°, as expected by symmetry; it is negative (downwards) 
for angles between 0° and 90°, with a maximal value of 
3 mN at an angle of about 40°. 

These measurements are independent of the initial ori- 
entation, even in the region of quickest variation (insert 
of Fig. |3J. This suggests that the results, obtained in 
transient regimes, would be the same if we could fix 9 to 
perform steady flow measurements. In fact, at a lower 
flow rate (25 ml min -1 ), we observe very similar tenden- 
cies, although more noisy (data not shown). It is thus 
natural to neglect the obstacle's inertia, and assume that 
the torque exerted by the flowing foam is exactly bal- 
anced by a friction torque (arising mainly from viscous 
dissipation in the capillary bridge between the ellipse and 
the top plate). Furthermore, the angular velocity \9\ is 
lower or comparable to 1° s -1 (Fig. hence the as- 

sociated Reynolds number a 2 \9\/v does not exceed 10. 
We can thus assume that the friction torque is propor- 
tional to the angular velocity 9, then Fig. 0Jd represents 
(up to an unknown multiplicative constant characteriz- 
ing the dissipation) the torque exerted by the foam. It 
displays a peak around 26°, compatible with Fig. [3] 

Fig. shows that the torque is negative for all pos- 
itive 9 values. Thus the only stable orientation of the 
ellipse is the parallel one, 9 — 0°. This contrasts with the 
Newtonian case, where long objects settle broadside-on 
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FIG. 3: Angle (x), drag (o) and lift (□) of the ellipse, versus 
time, for an initial angle of 76°. Insert: zoom on the region of 
quick variation of the angle. The data for three different initial 
orientations (x: 76°, +: 64°, A: 48°) are superimposed, by 
translating the time axis. The solid line is a fit to all data with 
a hyperbolic tangent profile, indicating that the maximum 
angular velocity is —2.1 ±0.1° s -1 for an angle of 26.4 ±0.1°. 



[13j . Note that in the case of a Stokes flow (without iner- 
tia nor elasticity), every orientation of the ellipse would 
be neutrally stable in an unbounded fluid 0], but that 
in the presence of boundaries, the parallel orientation 
is more stable ^f§. On the other hand, this is coher- 
ent with studies in other non-Newtonian fluids, where 
ellipses settle broadside-along under gravity in Oldroyd- 
B fluids or spherical particles aggregate vertically 
during sedimentation in shear-thinning fluids flB j. Ac- 
tually, the stable orientation of long objects under flow 
is determined by a competition between inertia and vis- 
coelasticity 0|> which have opposite effects. 

Fig. 0Jl shows that the lift is oriented downwards, as for 
a cambered airfoil, probably due to the positivity of the 
first normal stress difference 0| : this should therefore 
be valid for every viscoelastic fluid [TsT ]. 

It is worth noting that lift (Fig. EJt) and torque (Fig. 
0J}) are not maximal at the angle of 45°. This contrasts 
with the existing prediction of the torque exerted on an 
ellipse by a second-order fluid in potential flow , which 
predicts an angular dependence of the form cos 9 sin 9. 
We suggest two possible explanations for this discrep- 
ancy. First, the flow of foam is not potential, and even 
breaks the x — > —x symmetry between upstream and 
downstream |2l|. Second, second-order fluids might not 
be good models for foams, because they do not include 
yield stress. 

In yield stress fluids, viscoplastic models predict that 
the drag on circular obstacles is proportional to the ra- 
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dius of the obstacle 22] as long as the yield stress remains 
the preponderant contribution to the total stress. This 
agrees with experiments on circles @ , and this is compat- 
ible with the proportionality of the drag with the leading 
length of the ellipse (Fig. [2J . However, this scaling with 
the leading length docs not hold for a square obstacle, 
which experiences a drag independent of its orientation 
|f| for reasons we do not understand yet. In addition, 
any orientation of a square obstacle is neutrally stable 
in a flowing foam 0, whereas it would align its diago- 
nals streamwise and spanwise in a viscoelastic liquid, as 
reported in (lj) . 

To summarise, we are not aware of a single macro- 
scopic, continuous (viscoelastic or viscoplastic) model 
which can explain the whole set of experimental data. 

We now propose an elastic, microscopic model, to catch 
the main qualitative features of drag, lift and torque. 
We estimate the contribution F of the soap film tension 
(which determines the normal tensile elastic stress 0,0]) 
to the force on the ellipse. Since the foam is quasi-2D, 
each film separating two bubbles in contact with the ob- 
stacle exerts on it a force directed along the film; its 
magnitude is the line tension, A, which is twice the air- 
water surface tension 7, multiplied by the foam height 
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FIG. 4: (a) Drag (positive values) and lift (negative val- 
ues), (b) angular velocity 8 versus 8. Data correspond to the 
time-dependent experiments with four initial orientations: 76 
(same data as in Fig. |3J, 64, 48 and 18°. Also shown for 
lift and drag are data from the steady orientations, 6 — 0° 
and 90°. In (a), the plain curve is the fit by the spanwise 
dimension Q, and the dashed one is the fit by Eq. for 
the lift; both fits are up to a free prefactor. In (b), data for 
15° < 8 < 40° are noisy but compatible with the insert of Fig. 
|S| (open circle). The curve is the fit by Eq. (|5J, in arbitrary 
units. 



h, and a prefactor accounting for 3D geometry [l^. If 
the flow is quasistatic, the film is along the normal n to 
the surface of the ellipse (see Fig. \\$. The total force is 
thus a sum performed over the films in contact with the 
ellipse: F = Yl A/i. We do not model the contribution of 
the bubbles pressure, which is of the same order of mag- 
nitude, and is roughly proportional, to the contribution 
of the film tension |'2ll |23j • We do not model either the 
velocity-dependent forces and torque, originating from 
the viscous friction within the lubrication films between 
the ellipse and the surrounding bubbles. 

If the ellipse is much larger than the bubbles, we con- 
sider the distance between consecutive films along the 
ellipse as the continuous function /(a), a being the an- 
gle in the ellipse's parametric equation: X(a) = a cos a, 
Y(a) = fesina, and write the force F and torque C as 
integrals: 
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We then deduce the drag and lift as F x — Fx cos — 
Fy sin 6, and F y — Fx sin 9 + Fy cos 6, respectively. 

We must now model the function /, or equivalently, the 
deformation of bubbles around the obstacle. As already 
mentioned in 0] , this is strongly correlated to the local 
structure of the flow: if it converges towards the obstacle 
(leading side), it squashes the bubbles in contact, and / 
is high. Conversely, if the flow diverges from the obstacle 
(trailing side), it stretches the bubbles in contact, and / 
is low. Experimental images support this argument (Fig. 
0}, and, more precisely, lead us to set a phenomenological 
expression for /. Fig. ^ shows that the bubbles remain 
squashed over the whole leading side (/3 < a < n + /3 
with — arctan(6cot#/a) from elementary geometry); 
we thus assume that / takes a maximum value, /m over 
this interval. At the trailing side, Fig. ^ shows that 
the bubbles are progressively stretched up to a maximum 
close to the y — point, which is close to the angle 
a = —9 for simplicity. To reproduce this observation, 
we assume a piecewise affine variation of / from /m to 
a minimal value f m in the ranges —0<a<(3, and 
j3 — 7T < a < —6. The analysis of several images of the 
bubbles along the obstacle yields the following estimates: 
f m = 3.3 mm, and /m = 4.9 mm. Given the aspect ratio 
a/b = 1.6, we can calculate the drag, lift and torque from 
Eq. © (Fig. E|). 

For the drag, it turns out that the result from Eq. © 
is indiscernible (with 1% precision, up to a free prefac- 
tor) with Eq. Q; the agreement with the experimental 
data is thus quite good (Fig. QJi,). For the lift, we pre- 
dict the sign, i.e. explains the downwards lift: the tensile 
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stress is larger at the trailing edge where it contributes 
in average downwards (and downstream) for angles be- 
tween 0° and 90°, than at the leading edge, where it 
contributes upwards (and upstream). This confirms that 
the lift is dominated by the elasticity, as is the case for 
an airfoil [18|. Moreover, we predict correctly the an- 
gular dependence of the lift, and a maximum at angle 
40°, which agrees quantitatively with the experiments. 
For the torque, the agreement is qualitative: we predict 
its sign, the existence of a maximum at an angle smaller 
than 45°, and the stability (instability) of the parallel 
(perpendicular) orientations. 

The present model relies mainly on the coupling be- 
tween bubble deformation and flow. This argument has 
a very general validity: it explains the anti-inertial lift 
exerted by a flowing foam on an airfoil , and predicts 
quantitatively the drag on a circle on several decades of 
fluid fractions It also applies in 3D, as shown by the 
analogies between the 2D flow around a circle Q and the 
3D flow around a sphere H 0] . It is qualitatively insen- 
sitive to the presence of channel walls, both because this 
does not influence the convergence or divergence of flow 
close to the obstacle, and because of the very limited lat- 
eral extent of the influence of an obstacle for foams 0, Q 
compared to Newtonian fluids. 

Foams are often modelled as viscoplastic fluids, such as 
Bingham or Herschel-Bulkley models |24J. Such models 
describe yielding, which is occurring at the leading side, 
where the roughly constant amplitude of bubble defor- 
mation (Fig. ^) is a manifestation of yield strain. On 
the other hand, viscoelastic fluids such as the Oldroyd-B 
model often used for polymers |25j describe the delayed, 
elastic response of the bubbles, apparent at the trailing 
side through the progressive stretching of the bubbles 
(Fig. DO . 

Our phenomenological model captures both the cou- 
pling between strain and flow (with delayed response) 
and the saturation of deformation (yielding). It yields a 
good agreement with experimental data (Fig. 0}. Still, 
we can suggest improvements in three directions. First, 
the law assumed for / is a phenomenological description 
of observations. The next step would consist in predict- 
ing this function. This would require to quantify accu- 
rately the evolution of strain due to advection and plas- 
tic rearrangements of bubbles, which is predicted by re- 
cent models in the simple case of shear flow here, 
a generalization to more complex flows and geometries 
is required. Second, we could extend this model to de- 
scribe the velocity-dependent contribution to drag, lift 
and torque. This requires to quantify precisely the influ- 
ence of friction on the interaction between bubbles and 
obstacle boundaries |3j. Third, it would be useful to in- 
clude the effect of sharp angles, in order to understand 
why the drag on a square does not depend on its orien- 



tation. 
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